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Abstract 

Ultrafast laser excitation can induce fast increases of the electronic subsystem temperature. The subsequent electronic evolutions 
in terms of band structure and energy distribution can determine the change of several thermodynamic properties, including one 
essential for energy deposition; the electronic heat capacity. Using density functional calculations performed at finite electronic 
temperatures, the electronic heat capacities dependent on electronic temperatures are obtained for a series of metals, including free 
electron like, transition and noble metals. The effect of exchange and correlation functionals and the presence of semicore electrons 
on electronic heat capacities are first evaluated and found to be negligible in most cases. Then, we tested the validity of the free 
electron approaches, varying the number of free electrons per atom. This shows that only simple metals can be correctly fitted with 
these approaches. For transition metals, the presence of localized d electrons produces a strong deviation toward high energies of 
the electronic heat capacities, implying that more energy is needed to thermally excite them, compared to free sp electrons. This is 
attributed to collective excitation effects strengthened by a change of the electronic screening at high temperature. 

Keywords: Ultrashort laser; Femtosecond laser-excited electrons; First-Principles calculations; Electronic Fleat Capacities. 


1. Introduction 

Material response to intense laser excitation is the subject of important research activities, and recent advances 
have revealed the determinant role of primary excitation events. Their accurate comprehension is necessary to cor¬ 
rectly describe ultrafast structural dynamics lUBl^ phase transitions |13,01> nanostruture formation 0, ablation dy¬ 
namics HU , or strong shock propagation |@] . The interplay between ultrafast excitation and resulting excited material 
response still requires a comprehensive theoretical description for highly excited solid materials. Ultrashort laser irra¬ 
diation produces a strong inhomogeneous heating of electrons that rapidly exchange energy through electron-electron 
collisions, leading the electronic subsystem to a thermalized state which is intimately related to the electronic structure 
of the material. This fast heating of the electronic subsystem leads in turn to a significant electron-phonon nonequilib¬ 
rium, as the energy of the laser pulse is deposited before relaxation takes on and the material starts dissipating energy 
by thermal or mechanical ways. 

During intense laser exposition, the irradiated material undergoes successive stages of extreme constraints, start¬ 
ing from the inhomogeneous excitation of electrons, the swift rise of electronic temperature and electronic pressure 
after thermalization of the electronic subsystem on the tens of femtosecond timescale |@]. These rapid processes 
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are followed by the rise of the ionic temperature associated to thermally-triggered phase transitions, before ablat¬ 
ing or returning to ambient conditions from further energy dissipations occurring at larger timescale. The accurate 
knowledge of these rapidly evolving conditions, in intensity and time, is of importance as they modify material 
properties. Evolution of electronic and ionic temperatures, and T,- respectively, are generally calculated with the 
two-temperature model MM which involved transient parameters such as electronic heat capacity, electron-phonon 
coupling strength, thermal conductivity, electronic pressure and material absorptivity. 

The basic underlying assumption in the two-temperature model is that the electron subsystem consists of a ther- 
malized population of a certain amount of electrons, a condition easily fulfilled in ablation processes conditioned by 
high electronic temperatures ||3l. The energy evolution is determined by photon absorption, energy accumulation and 
thermal/mechanical transport (involving electron-electron and electron-phonon interaction) based on nonequilibrium 
dynamics. The effective number of charge carriers taking part in absorption, storage and dissipation of the energy 
has to be carefully evaluated to give a correct description of the evolution of the corresponding transient properties. 
A certain amount of free electrons responds collectively to the laser field through intraband transitions, associated 
to one-electron excitation corresponding to interband transitions. Modelling this absorption stage is rather complex, 
particularly for ultrashort laser irradiation as it initiates a strong charge disorder. The consequence of a nonthermal 
electron-ion population on the electronic energy storage is also of major importance. Finally, thermal and mechan¬ 
ical transports are often described by the concept of ’’free carriers”, since mobility is required for these processes. 
Electronic pressure has been seen to follow a free electron like behavior at high T^, even at solid density lll2n . The 
electronic heat capacity Ce plays an important role since it connects the quantity of absorbed energy to the rise of 
the electronic temperature of the system, regardless to energy dissipation processes. Electronic heat capacities have 
been derived from post-treatment of first-principles calculations where electronic excitation was not fully taken into 
account EH. Further, these quantities have been refined from density functional calculations dependent on the 
electronic temperature ifTil [fsll . 

In this paper, the electronic heat capacity of a series of transition metals is discussed in detail, first according to 
various modelling conditions, considering the effect of the exchange and correlation energy functionals and then the 
effect of semicore electrons. Ce are then discussed in the framework of a free electron approach in order to evaluate 
the impact of d electrons of transition metals on this crucial thermodynamic quantity. 


2. Calculation details 


The modelling of Al, Ni, Cu, Au, Ti and W metals is performed with the code Abinit lllbtl . which is based 
on p lane-waves description of the wavefunctions. Calculations are carried out within the density functional theory 
I ITIIiSD extended to finite electronic temperatures lll9n . Projector augmented-waves atomic data ll20l] are used to take 
into account the effects of nuclei and core electrons. A cutoff energy of 40 Ha is applied to restrict the number of 
plane-waves and the Brillouin zone is meshed with a30x30x30 k-point grid using the Monkhorst-Pack method 02111 . 
The local density approximation (EDA) and the generalized gradient approximations (GGA) are used with 
or without semicore electrons depending on the metals, in order to evaluate effects of semicore electronic states on 
computed properties. Computations are realized at the theoretical equilibrium lattice parameter of the crystal phase 
obtained at T, = 0 K. Thus, electronic temperatures from 10“^ to 10^ K are applied through a Fermi-Dirac distribution 
of electrons within a cold lattice. Such high Te are conjectural and interrogate about phase stabilities. As we need to 
achieve a high temperature asymptotic behavior in this study, we will assume here that time and spatial conditions are 
not fulfilled for phase transitions to occur. 

The calculated electronic heat capacity is derived from the variation of the internal energy E with respect to the 
electronic temperature as Ce = dEldTe- This thermodynamic quantity is computed for a series of metals at various Tg 
using EDA or GGA functionals with or without semicore electrons. Results are provided in Fig. [T]alongside the values 
of Lin et al. obtained from electronic structures calculated at Eg = 0 K Q. At a first glance, the agreement between 
their results and ours is rather good for the low temperature range, the deviations appear at Tg above 4 x 10"^ K, where 
electronic structures start to react significantly to the heating of the electronic subsystem. In particular, transition 
and noble metals having a c/-block within their valence band are sensitive to Fermi smearing, impacting d electron 
population, that induces a change of the electronic screening which in turn produces a significant energy shift of the 
(f-block HE. At high temperature, this produces a significant effect on the band structure, with consequences on 
the electronic heat capacities. The effect of exchange and correlation functional appears to be weak on computed Cg. 
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The LDA method generally provides lower equilibrium parameters than the GGA, leading to slightly higher electronic 
electronic heat capacities per unit of volume. 

The effect of semicore electronic states are also evaluated and Fig. [Hf) shows clearly that 4/ semicore electrons 
of W have an impact at high Tg. Despite the fact that the 14 / electrons lie on deep states, around 19 eV below 
the valence band, they are significantly impacted by the increase of Tg starting from 4 x 10“^ K. For Al, Ni and Cu, 
the highest semicore electronic states correspond to 2p ones, they are located 50 eV below the valence band and are 
weakly impacted, even at high Tg. Nevertheless, a small effect of semicore electronic states is noticeable, with stronger 
Ce values for Al and Ni and lower values for Cu, compared to calculations where semicore states are not considered. 
With first semicore electronic states located at -27 eV for Ti, a significant effect of semicore electronic states is likely 
effective at high Tg. A deeper investigation of the evolution of the electronic semicore states when Tg increases shows 
that they are as sensitive as the d band to changes of the electronic screening. Depopulation or population of the d 
band with Tg, induces either a decrease or an increase of the electronic screening which is signalled by changes of the 
Hartree energy, and which is effective on localized electrons of the system, from d electrons to semicore electrons. 
As an example, at 10^ K the energy of the electronic semicore states are shifted for -3 eV for 2p^ electrons of Ni, 
-8 eV for 2p® electrons of Cu, 2 eV for 2p® and -4 eV for 4/*^ electrons of W. Semicore states are not shifted for 
Al, since this metal does not have a d block in its valence band. The shifts of these deep electronic states follow 
the changes of the electronic screening as discussed in Ref. lll2n . This indicates that the evolution of the electronic 
screening, induced by changes of the electronic population of the d block, propagate to the localized charge density 
corresponding to semicore electronic states. 

Nonetheless, due to the low energy of electronic semicore states and their corresponding electronic densities more 
localized around the nucleus, semicore electrons thermalize slower than valence electrons 13 ■ According to this 
consideration, the following part of this study relies on PAW atomic data without semicore electrons, at least when 
they are not needed for a good description of general properties. This implies that the GGA functional is used for Al, 
Ni, Cu and the LDA for Au, neglecting the effects of semicore electrons. For Ti and W, the GGA functional is used 
with semicore electrons but without the 4/*^ electrons of W. 


3. Electrons involved in C. 


In the framework of a free electron model, Cg has a linear behavior at low temperature that can be modelled 
by a low temperature limit Cg - tP-N fUik^glgllef and saturates at high temperature to a non-degenerate limit - 
3l2NfnikB. Here, Nf represents the number of free electrons per atom, n, the ionic density, kp the Boltzmann constant 
and Sf the Fermi level. In Fig. |2 Cg is plotted for all the studied metals, and free electron like limits are presented at 
low and high temperatures. These limits are fitted to Cg with an effective number of electrons , provided in the 
figure, such as Nf - N'i/^ . At low temperature, the fit is performed up to an electronic temperature of 5 x 10^ K, while 
at higher temperature, the corresponding number of electrons is a rough estimation since the asymptotic limit is not 
reached for most of the studied systems. 

In order to facilitate the comparison between the number of effective electrons required to fit Cg and classical 


free electron numbers, their respective values are reported in Table [T] Nf from Ref. Il25ll are not designed to evolve 
with the increase of Tg due to the rigid definition of this number, which is based on the most extended orbitals 
of the electronic configuration of the isolated atoms. On the contrary, temperature dependent Nf{Tg) are obtained 
from extended-orbital overlapping considerations based on temperature-dependent electronic structure calculations 
II12I] . that allows electronic transfer between bands and Fermi-Dirac redistribution of electrons. Unsurprisingly, both 
free electron numbers remain quite close at low temperature, with a significant increase of Nf(Tg) compared to Nf 
when the temperature increases. This originates in the redistribution of electrons from localized d electronic states to 
delocalized electronic states of higher energies. 

As expected, both low and high temperature approaches perfectly catch the limit behaviors of the free electron 
metal Al, with an effective number of electrons equals to 3, in agreement with the numbers of free electrons expected 
for this element. However, the number of effective electrons, needed to correctly fit the complex evolution of electron 
heat capacities of transition metals, can be very high and somewhat unrealistic. At low electronic temperature, an 
important gap is observed between and expected Nf for transition metals. For W, the number of effective 
electrons equals the number of electrons in the valence band, while for Ti is even higher than available electrons 
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in the valence band. Reducing the temperature interval where the fit is applied does not modify significantly , 
except in the case of Ni where this number reaches 19 electrons to catch the particularly strong slope of Ce at - 300 
K ifl^ . In all cases, remains remarkably high compared to the number of free electrons classically used. 

Table 1. Typical free electron numbers Nf from Ref. and temperature-dependent free electrons numbers Nf(Tc) from Ref. fl^ . with estimated 
eff 

number of effective electrons Ng ' required to fit Ce in the low and high temperature regimes. 



Al 

Ni 

T. 

Cu Au 
, = 5x10^ 

Ti 

K 

W 

Nf 

3 

2 

1 

1 

2 

2 

Nf(Tg) 

3.0 

1.5 

1.9 

2.4 

1.3 

2.1 

n/^ 

3 

5 

II 

5 

10^ K 

6 

6 

Nf 

3 

2 

1 

1 

2 

2 

NfiTg) 

3.0 

2.9 

3.3 

4.2 

2.2 

3.5 

Nf 

3 

5 

7 

8 

6 

6 


From the low temperature to the high temperature regime, values evolve differently from a material to an 
other, remaining constant for Ni, Ti and W, while increasing for Cu and Au. These values still strongly exceed typical 
values of Nf, even if the temperature dependence is considered, signalling a collective effect of c/-electrons. 

4. Discussion 

The significant differences between and the typical values of Nf at low and high temperatures, demonstrate 
the inability of asymptotic limits derived from free electron approaches to grasp the complex behavior of transition 
metals, which relies on the presence of localized d electrons. Indeed, these localized electrons also contribute to the 
change of the electronic heat capacity, through the change of free energy in the system as Tg increases. The electronic 
heat capacity informs about the amount of energy needed to rise the electronic temperature of a given system. The 
fast deviation toward high values of Cg observed for transition metals, which is illustrated by values of n/-^ much 
higher than the classical values of Nf, implies that the energy required to heat the localized d electrons is higher than 
the energy needed to heat free electrons. Moreover, non-linear evolution of n/^ with the increase of Tg also indicates 
a complex impact of d electron contribution to Cg. Finally, previous studies showed that the change of d electron 
numbers with Tg modifies the electronic screening that in turn induces an energy shift of the r/-block affecting all 
localized d electrons HQ. Consequently, the effect of d electrons on Cg has to be considered as a collective effect 
of the overall localized d electrons. Thus, the number of electrons to be considered in the electronic heat capacity 
limits at low and high temperature should be reevaluated, according to an activity coefficient applied to the different 
types of electrons: 


Ng — (XfNf -b (XglNgl -b (XscNsc (1) 

where Nf represents the number of free electrons as defined before, Nj - A,, - Nf is the number of d electrons, with 
A,, being the number of valence electrons, and Nsc is the number of semicore electrons, a, stands for the activity 
coefficient for each category of electrons, af is set to the value of 1 for all metals as it corresponds to the free 
electron approaches. Since PAW atomic data used here exclude most of the semicore electrons and Nsc is set to 0. 
Considering that n/‘^ is equal to the values provided in Fig. I^it is then possible to deduce the activity coefficients 
of d electrons for all the transition metals at low and high temperatures from ad = (A^^ - Nf)INd- This is solved 
within a temperature dependence of free electron numbers and d electrons numbers that are presented in Table |2] The 
corresponding activity coefficients of d electrons are provided in Fig. |3] 

As expected, Od for Al remains equal to 0 due to the absence of d electrons. For Ni, Cu and Au, the coefficient 
values are in the range of 0.2-0.5, with an increase for Cu and Au and a decrease for Ni when Tg increases. These 
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Table 2. Temperature-dependent free electron numbers Nf{Te) from Ref. Cl and corresponding number of localized d electrons Nd{Te) = 
Nv — Nf{Te) in the low and high temperature regime. 



Tg = 5 X 10^ K 

Tg = 10^ 

K 


^f{Te) 

f^dCTe) 

Nf(Tg) 

NdiTg) 

Al 

3.0 

0.0 

3.0 

0.0 

Ni 

1.5 

8.5 

2.9 

7.1 

Cu 

1.9 

9.1 

3.3 

7.7 

Au 

2.4 

8.6 

4.2 

6.8 

Ti 

1.3 

2.7 

2.2 

1.8 

W 

2.1 

3.9 

3.5 

2.5 


significant activity coefficients of d electrons - that applies on all localized d electrons - imply that the response of 
these electrons is strong, even when the Fermi smearing weakly affects them at low temperature, as for Cu and Au. 
On the contrary, aj reach high values for Ti and W, around 2 and 1 respectively and slowly evolve with Tg. The 
strong values obtained for Ti and W originate from the fact that the Fermi level lies within the d block, even when Tg 
increases. Excitation effects directly involves population or depopulation of localized d electronic states with strong 
impact on Cg. With a Fermi energy located on the top of the d block at low electronic temperature, a similar effect 
should also be observed for Ni, but this is actually hidden by the fit on large range of electronic temperatures. As we 
will see later, this phenomenon becomes visible when the fit is performed at lower temperatures. For W, once aj is 
obtained from calculations that neglect the effect of semicore electrons, it is then possible to use this value to derive 
Qsc from the Cg curve obtained by taking into account the 4/'^ electrons (GGAsc curve in Fig. [Ilf))- This leads to 
a new fit in the high temperature limit, with being equal to 9. Accordingly, Osc equals to 0 at low temperature 
and reaches the value of 0.2 at high temperature, confirming the effect of these electrons on Cg, assuming they are 
therm alized. 

Since Cg is supposed to evolve linearly in the low temperature regime, it is often characterized by the electronic 
heat capacity coefficient y - CgITg. Values at 5 x 10^ K are provided in Table[3] where a relative good agreement is 
found between calculated values and the low temperature free electron limit y' based on . However, for Ni the 
linear evolution of Cg is not respected, and a strong change of the slope occurs at low temperature. At Tg equals to 300 
and 10^ K, the corresponding y values reach 1183.3 and 631.1 J m“^ respectively. In the low temperature limit, 
this involves very high values of , 19 and 10 respectively, leading to activity coefficients of d electrons of 2.1 and 
1.0 respectively. This strong change of Cg slope at low temperature corresponds to a major difference in the behavior 
of Ni compared to Cu and Au. This difference of behavior originates from the position of the Fermi level, at the top 
of the d block for Ni while it is higher for Cu and Au, due to difference in their number of valence electrons. Thus, 
at low temperature, depopulation of the d block is direct for Ni, while the d block is not impacted yet for Cu and Au. 
When Tg increases, the Fermi level shift toward energies higher than the d block for both elements, inducing a strong 
changes of in case of Ni, that progressively tends toward the values of Cu and Au. As in the case of Ti and W, a 
direct impact of the Fermi smearing on the d block has a strong effect on Cg which is signalled by important activity 
coefficient of localized d electrons. To a lower degree, a similar deviation is observed for Ti around Tg - 10^ K, with 
a y coefficient reaching the value of 466.3 J m~^. This is fitted by a value of equals to 8 and leads to an 
even higher ad with the value of 2.5. For these elements, the non-linear evolution of Cg at low electronic temperatures 
leads to strong modifications of subsequent derived quantities. Thus, beyond these general observations, the evolution 
of these activity coefficients remains difficult to be evaluated accurately and values cannot be extrapolated to other 
elements or temperature regime. 


Table 3. Coefficient of electronic heat capacity (y = CgITg, J K ^ m ^), at 5 X 10^ K. 



Al 

Ni 

Cu 

Au 

Ti 

W 

r 

110.4 

320.0 

243.9 

141.1 

340.9 

266.8 

r' 

96.2 

309.3 

262.5 

172.8 

328.7 

342.2 
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5. Conclusion 

In this paper, the electronic heat capacities are obtained from first-principles calculations performed at finite elec¬ 
tronic temperatures. The effect of exchange and correlation functional on this thermodynamic quantity is first dis¬ 
cussed. At high electronic temperature a small difference is observed between LDA and GGA functionals that was 
attributed to different equilibrium lattice parameters. The impact of semicore electrons, is also tested and demon¬ 
strated at high temperature for 4/^"^ electrons of tungsten. This is dependent on the energy depth of the corresponding 
semicore electronic states and is mostly negligible even at high temperatures, except for W. Assuming negligible ef¬ 
fects or lower thermalization times for these electrons, effect of semicore electrons were neglected in the present study. 
Electronic heat capacities are also compared to previous theoretical predictions obtained from electronic structures 
not relaxed with the electronic temperature, with a very good agreement up to 4 x 10“^ K. Beyond this temperature, 
the response of the electronic structure exhibits differences. 

The obtained electronic heat capacities are then discussed in the framework of low and high temperature limits of 
a free electron approach. Both are based on a free electron number, which is modified up to an effective number of 
electrons in order to fit the theoretical values of Cg- This approach applies correctly to A1 as a free electron metal, but 
a significant overestimation of free electron numbers is found for the transition metals. Assuming a collective contri¬ 
bution of localized electrons on the Cg evolution, d electrons where considered associated to an activity coefficient. 
The high values obtained for this coefficient tend to corroborate the fact that d electrons respond collectively to the 
electronic excitation. A strong effect is observed from the impact of Fermi smearing on d electrons, as showed for Ni 
at low temperature and Ti and W in the whole range of temperature. For W, at very high temperature, an important 
contribution of semicore electrons 4/ is also deduced from this approach, assuming that they are thermalized. At high 
temperature, the change of the electronic screening is found to affect significantly the localized d electrons, with a 
global shift of the ^f-block, strengthening the collective behavior of these electrons. 

Additionally, the non-linear effect of d electrons, which is shown by an important dispersion of their activity 
coefficient, illustrates their non-free behavior and their interactive nature. The large numbers of electrons needed to fit 
Ce evolution in the low and high temperature regime is a signature of localized electron influence, mostly d electrons 
but also semicore electrons at very high temperature. Consequently, an unique free electron number cannot encompass 
the complex evolution of the electronic heat capacity and localized electrons, especially d ones, have to be considered. 
This shows that the free electron number is only one of the components on the energy storage of an irradiated material, 
and that localized electrons contribute as well, to various degrees depending on the electronic temperature. 
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Figure 1. (Color online) Evolution of the electronic heat capacity with T^, for Al (a), Ni (b), Cu (c). An (d), Ti (e), W (f), using 
different EDA and GGA exchange and correlation functionals. GGAsc and LDAsc indicate that the effect of semicore electrons 
are taken into account. The dotted black curves represent data of Lin et al. Q. 
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Figure 2. (Color online) Evolution of the electronic heat capacity with for A1 (a), Ni (b), Cu (c), An (d), Ti (e), W (f), and its 
corresponding low and high temperature limits. The number of electrons used to fit the theoretical curve in both models is provided. 



Figure 3. (Color online) Histogram representation of activity coefficients of localized d electrons. For each material, the left 
histogram corresponds to the low temperature case while the right one (patterned) represents the high temperature case. 
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